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Abstract 

We reanalyse the LSND neutrino oscillation results in the framework of the Projected Quasiparticle 
Random Phase Approximation (PQRPA), which is the only RPA model that treats the Pauli Principle 
correctly, and accounts satisfactorily for great majority of the weak decay observables around ^^C. We 
have found that the employment of the PQRPA inclusive DIF ^^C(z^e) e~)^'^N cross-section, instead of 
the CRPA used by the LSND collaboration in the z/^ — > I'e oscillations study of the 1993 — 1995 data 
sample, leads to the following: 1) the oscillationprobability is increased from (0.26 it 0.10 it 0.05)% to 
(0.33 lb 0.10 lb 0.13)%, and 2) the previously found consistence between the (sin^ 26, Am^) confidence 
level regions for the z/^ — > i'e and the D(, oscillations is significantly diminished. These effects are 

not due to the difference in the uncertainty ranges for the neutrino-nucleus cross-section, but to the 
difference in the cross-sections themselves. 
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Several recent experiments 

mmm 

strongly suggest that neutrinos oscillate. This 
means that a neutrino of a certain flavor (e.g. z/^) transforms as it propagates into a neutrino 
of another flavor (e.g. z/g)? violating the conservation of the lepton number. For this to happen, 
the simplest and most widely accepted explanation is that neutrinos have masses and mixing. 
There are evidences of transitions for three different Am^: ~ 8 x 10~^ eV^ (solar), ~ 2.5 x 10^^ 
eV^ (atmospheric) and ~ 1 eV^ (LNSD), which cannot all be understood in the context of three 
neutrino oscillations. Normally, the Liquid Scintillator Neutrino Detector (LSND) results are 
not taken into account when fitting neutrino oscillation data. Nevertheless, one has to try to 
understand the real significance of the LSND measurements, specially because Am^ ~ 1 eV^ 
is of particular interest to astrophysics and cosmology. The LSND experiment took place over 
six calendar years finding evidence for the appearance of electron-antineutrinos Ue at the 3.3a 
level ^ 0], and at lesser significance they have observed as well hints for the appearance of 
electron-neutrinos i/e 

. The first of these signals is the main LSND result and the second 
weaker signal was used as a consistency check. The positive results in both channels were 
interpreted in a two-flavor framework as transitions between the weak eigenstates (z/^) and 
{^e) driven by masses and mixing. In fact, quantum mechanics dictates that in this case the 
normally observed weak eigenstates (z^^, Ve) can oscillate between each other with probability 

P.^_., = sin\2e) sin^ (l.27 Am^^) , (1) 

if they are composed of a mixture of mass eigenstates (z/i,z^2)- Here 9 is the mixing angle 
between the mass and flavor bases, Am^ = m\ — m\ is the vi and V2 mass squared differences 
in eV^, Ly is the baseline, the distance in meters travelled by the neutrino from the source to 
the detector, and Ey is the neutrino energy in MeV. 

The combination of the LSND data with other compelling evidences for neutrino oscilla- 
tions, stemming from atmospheric ^|, solar ^, KamLAND ^ reactor, and K2K [?| accelerator 
neutrino experiments, cannot be adequately explained in the standard three-neutrino picture 
with CPT conservation, and this issue is considered to be a big challenge for neutrino phe- 
nomenology y, y, lioi. Models with four light neutrinos (the extra neutrino being sterile) fll| 
or CPT violation (1^ with three neutrinos have been proposed to accommodate all neutrino 
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data. However, in both cases, recent analyses show that neither scenario provides a satisfactory 
description of the data 

HQ. 

In the LSND experiment the neutrinos come from the decay of vr^ in flight {decay in 
flight, DIF), whereas the neutrinos z/g and the antineutrinos z/^ come from the decay of fi^ at 
rest {decay at rest, DAR), i.e., 

i 

e+ + Z/e + u^. 
DIF DAR 

The search for the DAR z/^ — > z/g oscillations involves the measurement of the reaction 
p{i'e, e~^)n, which has a large and well known cross section. The events are identified by detecting 
both the e"*" and the 2.2 MeV 7-ray from the reaction p{n, '~f)d. On the other hand, the signature 
for the DIF z/^ v^, oscillations is marked by the presence of an isolated high energy 

electron (60 < E^^^ < 200 MeV) in the detector. It is produced by the charge-exchange 
reaction ^^C(z/e, e~)^^A^, which takes place throughout the tank, the cross section of which 
is, as yet, not well established. The lower and upper energy cuts for E^^^ were chosen in such 
a way as to be above the Michel electron end point of 52.8 MeV and below the point where the 
beam-off background starts to increase rapidly and the signal becomes negligible. 

There are two LSND studies of the DIF — > z/e oscillations. The first analysis was done 
on the 1993 - 1995 data sample [2], which gave a total of A^°f = 18.1 ± 6.6 ± 4.0 oscillation 
events, corresponding to a transition probability 

P;;^,^ = (2.6±1.0±0.5)xl0-^ (2) 

when the cross-section cxe predicted by Kolbe et al. within the Continuum Random Phase 
Approximation (CRPA) is used jl5|. In the second search, the 1996 — 1998 data sample |^ 
was included as well, with reduced DIF fiux and higher beam-off background compared to the 
1993 — 1995 data. The reason for this modification lies in the fact that in this study first priority 
was given to the DAR z/^ Vf. oscillations, which have been analysed jointly. Moreover, for 
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has been employed in this occasion two different models. Namely, the shell model (SM) 
estimate, done by Hayes and Towner for the DAR region, and a relativistic Fermi gas 

model for the DIF region. The resulting total excess was N°^'^ = 8.1±12.2±1.7 events, yielding 

P;;^,^ = (1.0±1.6±0.4)xl0-l (3) 

The aim of the present work is to explore the role played by these nuclear structure effects in 
the delimitation of the neutrino parameters for the DIF oscillations.^ 

Before proceeding, and to make more clear the objective of the present work, it is convenient 
to discuss briefly the fiux-averaged exclusive cross sections 

= / dE,a,{E, = E,-/\M)^,{E,) (4) 

•J 

and the inclusive cross sections 

ar = / dEMEu)MEu), (5) 
Jo 

where 

ae{E,) =J2ae{Ee = E,-ujj,Jjy, i = e,fi. (6) 







The spin and parity dependent cross section a£{Ei, JJ) is given by [17|, (2.19)], uj^, are the 
excitation energies in ^^N relative to the ground state in ^^C, and A = uj-^^+ = 17.3 MeV. The 
energy integration for electrons is carried out in the DAR interval mg + cuj^ < A < E^f'^ = 
52.8 MeV, and for muons in the DIF interval up to + cuj^ < A^f < E^''^ = 300 MeV. ^ 
^e{Ey) is the normalized neutrino flux; for z/g it is approximated by the Michel spectrum, and 



for z/^ that from Ref. UM was used. 



The experimental data for the exclusive and inclusive cross sections, given in Table lU show 
that the DAR and DIF processes are of quite different nature: while the first one is dominated 



^ Accurate knowledge of the i' cross-section, and the related observables, plays an important role for the next 
generation of experiments. Various target nuclei, like C, O, Fe, Ar, Pb, • • •, are normally (and presumably 
will be) employed to provide the detector mass. 

^ In order to invert the summation on J J and the integration on dE,y, we have extended the lower limit of 
integration in Q from mi + luj^ to zero by defining (7i{Ee = — ^Jj, J]) = ^ f'^'' ^ "^^^ 
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TABLE I: Calculated and experimental flux-averaged exclusive o^^, and inclusive ajf^ cross section 
for the ^'^C{i'e,e~y^N DAR reaction (in units oflO""^^ cm^) and for the ^'^C{u^, n^^^N DIF reaction 
(in units of 10"^'^ cm^). The CRPA calculations |l5| were used in the flrst LSND analysis on the 1993- 
1995 data sample 1^, and the SM calculations from Ref. in the second LSND oscillation search 



The listed PQRPA results correspond to the calculations performed with the relativistic corrections 
included One alternative SM result as well as the RPA and QRPA results from Ref. jl^ are also 
shown. 
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8.1 


18.6 


0.59 


13.0 


SM [19] 


8.4 


16.4 


0.70 


21.1 


RPA 


49.5 


55.1 


2.09 


19.2 


QRPA [19] 


42.9 


52.0 


1.97 


20.3 


Experiment 










Ref. [21]] 


9.1 ±0.4 ±0.9 


14.8 ± 0.7 ± 1.4 






Ref. [21J 






0.66 ±0.1 ±0.1 


12.4 ±0.3 ±1.8 


Ref. [22] 


8.9 ±0.3 ±0.9 


13.2 ±0.4 ±0.6 






Ref. [23] 






0.56 ±0.08 ±0.10 


10.6 ±0.3 ±1.8 



in proportion of 2/3 by the Gamow- Teller (GT) transition to the ground state 1^ in ^^A^, 
the second one populates almost entirely the excited states through the forbidden transitions. 
It is quite a difficult task for the nuclear structure models to describe both cross sections 
simultaneously. 
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The SM treats correctly the Pauh Principle within the p-shell, which is crucial for the correct 
distribution of the GT strength, whereas the predictions for high-lying states are less certain 
because of the truncation of the model space. In fact, the SM calculation performed by Hayes 
and Towner [3l reproduces fairly well several data. But, in a later SM study, Volpe et al. j|19|] 
noted that this concordance could be an artifact because the employed model space was not 
large enough to exhaust the charge-exchange sum rules. More, the same authors have shown 
that when a more extended space is employed the SM cross sections are increased exceeding 
the experimental LSND result. 

The RPA like models include high-lying one-particle one-hole excitations, but very frequently 
completely fail to account for the amount and distribution of the GT strength as can be seen 
from Table HI This is the reason why the CRPA is unable to explain the weak processes 
(/3-decays, /i-capture, and neutrino induced reactions) among the ground states of the triad 
ji2^ i2(^ i2jYj. ^ rescaling factor of the order of 4 is needed to bring the calculations and 
the data to agree [13], and a subsequent ad hoc inclusion of partial occupancy of the pi/2 



subshell reduces this factor to less than 2 



2J, |25|. It is still more relevant here that the 



CRPA overestimates the inclusive ^^C(i^,,,u )^^A^ cross-section with v.. coming from the DIP 

n n 

of 7r+ by about 50% |2l| or more [2^, because one can assume that the DIP ^^C(z/e,e )^^A^ 
cross section, which gauges the — >■ z/g oscillations, is affected in the same proportion. This 
assumption comes from the universality of the weak interaction and was done in the first LSND 
analysis ^ 

Thus, it might be interesting to reanalyse the LSND results in the framework of the Projected 
Quasiparticle Random Phase Approximation (PQRPA) ^|, which is the only RPA model that 
treats correctly the Pauli Principle, explaining in this way the distribution of the GT strength. 
To achieve this it was imperative both: a) to include the BCS correlations, and b) to perform 
the particle number projection. Under these conditions most of the weak decay observables 
around ^^C are within 20% of the PQRPA predictions. This happens, for instance, with: 1) 



^ Since the work of O'Connell, Donelly and Walecka |3| we know that electron and muon cross sections differ 
for low neutrino energy, but tend to merge for high neutrino energy. 
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the B(GT)-values to ^^A^ and ^'^B, 2) the exclusive inuon captures to the 1^, 2^, \i and li 
states, as well as the inclusive muon capture in ^^B, and 3) the exclusive cross sections o^^'^ 



and and the inclusive cross section a^^'^ |l7l. l27l|. The only exception is the inclusive cross 
section, a"'^* = a™'^ — o^^'^, for which the PQRPA value, 10.5 (in units of 10~^° cm^), is more 
than 100% larger that the experiment result, 4.3 ± 0.4 ± 0.6 j22]. From the nuclear structure 
point of view the theoretical evaluation of this quantity is a peculiarly delicate and subtle issue 
and therefore deserves a special comment. In fact from Table VI in Ref. [3] it can be seen that 
in the PQRPA case a™"^* is build up from the interplay of GT strength not contained in the li" 
state, the Fermi (F) transitions to the O"*" states, and the first forbidden transitions to the 1~ 
and 2^ states. All these quantities are relatively small and evaluating them precisely is a very 
difficult task. Then one should not be surprised by the most recent SM calculation which 
yields a result^ (ajf"^* = 8.3) which is twice as large as that obtained in the previous SM study: 
g.mc* ^ ^^Qj_ rjj^g CRPA result a'^^* = 6.3 llj, very likely does not contain any GT and F 
strengths as it should, and therefore, in this case, the agreement with the experiment could be 
accidental. 

We will limit our attention only to the 1993 — 1995 data sample P], which, as mentioned 
before, yields a more defined signal for the oscillation events. The experimental oscillation 
probability can be written as 

pexp ^ (^^^e) /yN 

where the z/g fiux (from now on) is defined as 

= < + (8) 



with the fiuxes and coming, respectively, from the DIF decays 7r+ ^ e"^ + z/g and 
/i+ e+ + z/e + v^,. fn = (9.23 x 10^^) ■ (5.4 x 10^°), with the first quantity being the number 
of protons on target (POT), while the second one is the fiducial volume (number of molecules 
of CH2 in the detector tank). = N"'" + N^s = 27.7 ± 6.9 is the total number of beam- 
excess events measured by LSND, and e is the event selection efficiency. The averaged inclusive 



7 



cross-sections are 

(a<l>,,) = J2 ^e{E, = E,- uj^, JJ)^,,dE,; i = e, fx, (9) 

where Efn = 60 MeV + ujj and Ej'j = 200 MeV + uj^ . In order to simplify the numerical 
calculations which follow, instead of using the exact equations (jH)), we will employ here the 
approximate ones: 

((T$,,) = r ae(^.)$.,rf^.; i = e,fi, (10) 
Je< 

where a^iE^) is given by ®, and E< = 60 MeV + A, and E> = 200 MeV + A. We have 
verified numerically that the equations (fTUj) reproduce the equations Q up to a few per cent. 
The neutrino fluxes ^^^,^1^ and were adopted from the Ref. [3]. The CRPA and PQRPA 
results for ae{Ey), (Te(-Ej/)$j/^ and ae{Ey)^ve are confronted in Fig. [TJ as a function of Ey. 

The systematic error associated with the PQRPA cross-section is taken to be 20%, based 
on our theoretical uncertainties (see Tables V, VI and VII in |l2i|), and agreement between 
measured data and theoretical predictions for the weak decay observables involving the ^^C 
nucleus 27 1. Therefore, considering the same uncertainties as in the LSND search in 



the selection of e (12%) and in the flux (15%), we end up with a total systematic error of 
28%, which yields N"'"" = 21.5 ± 6.6 ± 8.5. In this way the PQRPA result for the oscillation 
probability turns out to be: 

P;;^,, = (3.3±1.0±1.3)xl0-l (11) 

The difference when compared to the CRPA result Q is due to the difference in the electron 
cross-section, as evidenced in Fig. [TJ 

In order to determine a confidence region in the (sin^ 29, /S.m?) parameter space we proceed 
in the same manner as in Ref. Q]. First we rearrange the data for energy distribution of the 
excess events (see 0, Fig. 29]) in four equal energy bins Ny{i), as shown in Fig. |21 Next, we 
minimize the function 



:i2) 
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FIG. 1: (Color online) Comparison between the CRPA and PQRPA results for: (Te{Ey) in units of 
10"^'' cm^ (upper panel), and, in units of 10~^^ POT MeV for ae{E^)^iy^ (middle panel) and 
ae{E^)^^^ (lower panel). 



where N^{i) = N°'%i) + N^^ii), with 



J^OSCI 



= eU I a{E,)n{E,)^ll{E,)dE,, 



(13) 



where Pv^u^ is a function of E^, sin^(26') and Am^, and is defined in ((H). We include the 



resolution function 



7^(^.) 



27re(E, 



E> 



exp 



1 ( k-e; 

2 \ e{E,) ^ 



dE'. 



(14) 



which takes into account the probability for finding the electron inside the window of detection, 
with e{Ey) = O.OGE,^ being the experimental energy resolution 
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^14^ 

l3 i2[. 
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160 180 200 
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FIG. 2: (Color online) The energy distribution (in 4 energy bins) of the LNSD excess events, N^{i), 
together with the corresponding experimental errors 5Ny{i) (vertical lines) and the energy intervals 
Ei,{i) (horizontal lines). The theoretical CRPA and PQRPA values for the expected background 
events, N^^{i), are shown as well. 

I — I 

To set the confidence levels (CL) we used the raster scan method [29]: for each value of 
Am^, a best fit is found for svo?26. At each Am^, is calculated as a function of sin^ 2^. 
The ID confidence interval in sin^ 20 at Am^ is composed of all points having a within 3.84 
of the minimum value (3.84 is the two-sided 95 % CL for a distribution x^ with one degree of 
freedom). The confidence region in the (sin^ 26, /S.w?) is the union of all these intervals. 

Our (sin^ 29, Am^) oscillation parameter fits for the DIF channel u^, corresponding to 



both the CRPA 



1^1 and PQRPA 



13] cross-sections, are shown in Fig. IHl along with the favored 
regions for the LSND DAR measurement for t'^ — i^e In order to better understand the 
consequences of using different cross-sections, the confidence regions obtained with (lower panel) 
and without (upper panel) inclusion of the systematic uncertainties, are displayed separately. 
In the calculation with the CRPA cross-section these uncertainties are taken to be the same 
as in the LSND search P], i.e., of 22% for the positive side, which shifts the parameter space 
downwards, and of 45% for the negative side, which shifts the parameter space upwards. On 
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FIG. 3: (Color online) Regions in the neutrino oscillation parameter space. In the upper panel the 
results for — > oscillations without the inclusion of the systematic uncertainty are shown, while 
the lower panel shows those with the uncertainty included, as described in the text. 

the other hand, an uncertainty of 28% is used for both negative and positive side, when the 
PQRPA cross-section is employed. 

We see that, when the systematic uncertainties are considered, the CRPA 95% CL region 
fully comprises the 99% CL region for the — > Ve oscillations, which is in essence the result 
obtained by the LSND collaboration 0], Contrarily, this does not happens in the PQRPA case 
where the overlapping between the two regions is only marginal. It is important to stress that 
the v^j, — > z/g region is dragged towards the z/^ — > Ve region by the positive side uncertainty. 
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while the role played by the negative side uncertainty is of minor importance. For the sake 
of completeness the result of the joint z/^ —>■ v^. and z/^ Vf, oscillation parameter fit over 
(sin^ 2^, Am^) plane for the complete 1993 — 1998 data sample 3], is also displayed in Fig. El 

In summary, we have found that the employment of a smaller inclusive DIF ^^C(z/e, e~)^^A^ 
cross-section, than the one used by the LSND collaboration [3] in the z/^ Vf. oscillations study 
of the 1993-1995 data sample, leads to the following consequences: 1) the oscillation proba- 
bility Pl^^y^ is increased, and 2) the previously found consistence between the (sin^ 2^, Am^) 
confidence level regions for the Vy, Ve and the — > v^. oscillations is diminished. More, these 
effects are not due to the difference in the uncertainty ranges for the neutrino-nucleus cross- 
section, but to the difference in the cross-sections themselves, and are quite significant when 
the PQRPA is used instead of the CRPA. Thus, precise knowledge of the nuclear structure 
involved in the z/-nucleus cross- sect ion, could play an important role in the delimitation of the 
neutrino parameters for the DIF z/^ v^, oscillations. 
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